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ABSTRACT 

We have measured the proper motion (PM) and kinematics of the ansae in 
NGC 7009 using high dispersion echelle spectra and archive narrow band HST 
images. Assuming that the ansae are moving at equal and opposite velocities 
from the central star we obtain a system radial velocity of —53 ± 2 km.s~^, the 
eastern ansa approaching and the western ansa receding at Vr = 5.3 ± 1 km.s~^ 
with respect to this value. The PM of the eastern ansa is 28 ± 8 mas.yr"^, which 
with our weighted distance to NGC 7009 of 0.86 ±0.34 kpc gives Kxp = 112 ±32 
km.s~^. The electron temperature and density in both ansae were determined to 
be Te ~ 9000 ± 400 K and ~ 2300 ± 400 cm^^. The dynamic age of the ansae 
is ~ 925 ± 260 yrs. and the implied PM of the central star is fics = 1 ± 0.5 
mas.yr~^. This is in qualitative but not quantitative agreement with previous 
work. 

Subject headings: astrometry — ISM: kinematics and dynamics — planetary 
nebulae: individual (NGC7009) 



1. Introduction 

Most low and intermediate mass stars become White Dwarfs after a brief phase as 
Planetary Nebulae (PNe) caused by heavy mass loss on the AGB (Iben & Renzini 1983). 
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Symmetrical and mildly asymmetrical PNe are explained by versions of the interacting wind 
model (Kwok et al. 1978) but the mechanisms causing the strong asymmetries observed in 
the majority of PNe are still unexplained. That both binaries and magnetic fields probably 
play a role is now generally accepted. For an overview of formation mechanisms, models, 
and asymmetries see Kastner et al. (2000). 

The wealth of different PNc shapes is shown in Balick (1987) and Schwarz et al. (1992), 
and features have been described in detail by many authors: e.g. point symmetry (Schwarz 

1993) , bipolarity (e.g. Corradi & Schwarz 1995), multipolarity (Sahai & Trauger 1998), 
BRETS (Lopez 1997) and Fast, Low Ionization Emission Regions or FLIERS (Balick et al. 

1994) arc all used to indicate the various observed morphologies. 

NGC7009 has long been known to have a complex morphology. Aller (1941) described 
the outer features as ansae (handles in Latin), and work by Reay & Atherton (1985), and 
Balick et al. (1987) show that these ansae are expanding near the plane of the sky at high 
velocities (~ 10^ km/s) relative to the nucleus. 

Liller (1965) used photographic plate material to analyze the expansion of NGC7009, 
including the ansae, for which they obtained 16 ± 3 mas.yr"^. 

In this paper we present the kinematics and proper motion of the ansae using observa- 
tional material from the HST and CTIO 4m echelle spectrograph, and show that the previous 
results are qualitatively but not quantitatively correct. 

2. Observations and data reduction 

We extracted two sets of archive HST images with the "On the fly" option. The WFPC2 
sets were taken on 1996.04.28 (Bahck et al. 1998) and on 2001.05.11 (Palen et al. 2002) in 
an Ha/ [Nil] filter. 

The 2001 images are shifted 20 " east from the 1996 images. As the latter covered the 
entire nebula, in the former the west ansa disappears from the field of view. The two images 
showed in Figure 1 are an average (for cosmic ray rejection) of each set. 

We took CTIO 4m echelle spectra on 2002.07.29 with a spectral range between 410- 
720nm, and a slit of 1.2" x 6.6". The mean seeing was 1.4". Flux cahbration to 10% was 
done with 58 Aql (HR 7596, Hamuy et al. 1992). The spectra were reduced in the usual 
way, using IRAF, and resulted in an rms wavelength error of 0.00037nm, with residuals of 
< 0.00075nm. Velocity resolution was 3.7km.s~^.pix~^ or 10.2 km. s~^ for our slit width at 
650nm. 
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Fig. 1.— Average of the set of HST images taken in 1996 (left) and 2001 (right). North is 
down-left and east is down-right, and each image is 120" by 120". There is a shift of 20" 
in the east-west direction between them, which makes the western ansa disappear from the 
2001 images. The two background stars marked with circles were used together with the CS 
as reference points for proper motion calculation. 

3. Results 

3.1. Distance &; proper motion 

The distance to NGC 7009 is uncertain, but we computed the weighted average of the 
14 available (Acker et al. 1992) values to be ci = 0.86 ± 0.34 kpc. 

The measurement of the central star's position was made with the IRAF task imcentroid. 
For the measurement of the position of the ansa, we confined it inside a box of 5" square and 
computed the centroid, using the imcentroid task and measured their displacements with 
respect to the central star (CS) and two field stars (shown in Figure 1) in two ways: the 
displacements were measured in all the images, and an average displacement was obtained 
for each epoch; next, the displacements were computed from the averages of the images for 
each epoch. The centroid positions given in pixel values were converted to RA and DEC 
using the metric task, within the stsdas package. This task corrects for the known geometric 
distortions between different CCDs. The differences in the results are an order of magnitude 
smaller than the errors derived from centroid calculation and measurements between different 
chips. The various differential images are shown in Figure 2. 
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Fig. 2. — Difference between the combined 2001 and 1996 images, normalized by the average 
flux of the three reference points chosen: central star and the two field stars. The panels 
show differences wrt: central star (upper-right), field star 1 (lower-left), field star 2 (lower 
right) and the adopted least squares fit (upper-left). Each image is 43" by 24". 
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Tabic 1: The relative separation between central star and eastern ansa for each of the images. 
1996 images 



Dataset 


UT of exp. 


Relative separation 


(") 


u32e0308t 


20:48:16 


24.85 ± 0.04 




u32e0309t 


20:56:16 


24.84 ± 0.04 




u32e030at 


21:04:16 


24.86 ± 0.04 




Average 


20:56:16 


24.85 ± 0.06 




2001 images 








Dataset 


UT of exp. 


Relative separation 


(") 


u5hc6002r 


19:05:14 


25.01 ± 0.04 




u5hc6003r 


19:13:14 


25.00 ± 0.04 




u5hc6004r 


19:24:14 


25.00 ± 0.04 




u5hc6005r 


19:31:14 


25.02 ± 0.04 




Average 


19:18:29 


25.01 ± 0.05 





The displacements between the CS and the ansae are: 24.85±0.06"& 25.01± 0.05"(E; 
1996 & 2001) and 26.77± 0.06"(W; 1996 only). The detailed measurements with their 
respective errors are shown in Table 1. Note that the adjusted standard deviations for 1996 
and 2001 results are 0.01 and 9.6.10"^ respectively. 

Even though this is a good determination for the angular displacement of the eastern 
ansa, it is also possible that the central star is moving. Therefore, it is important to have 
other reference points. The circles in Figure 1 indicate the position of the two field stars 
present on the images that could be measured. The angular displacements of the ansae in 
each frame are Hsted in Table 2. These measurement are all equal within the errors. 

The value we adopted for the displacement was the least squares shift of all measure- 
ments: 0".14± 0".04. The time between images is 158883733s or 5yrs. computed from 
the start time of the exposures, yielding a proper motion of the E ansa of 0".028 ± 0".008 
yr~^, and the tangential velocity is Yexp — 112 ± 32 km.s~^ for d = 860pc. 

We note that the rim of the inner nebulae also shows differential motion in more or 
less the same direction as that of the ansae between the 1996 and 2001 images. We do not 
attempt to compute the proper motion as this rim is likely to be a moving ionization front 
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Table 2: Angular displacements of the east ansa as measured from the combined 1996 and 
2001 images using the reference points specified (Upper star and Lower star refer to Figure 1). 



Referonce Point Angular dis]:)lacement ('') 

Central star 0.17 ± 0.10 

Upper star 0.16 ±0.10 

Lower star 0.10 ±0.02 

Average shift 0.13 ± 0.04 

Least squares shift 0.14 ± 0.04 



with a very different velocity than that of the bulk motion of the gas (weak R-front, Mellema 
(2003)). This is unhkely to be the case for the ansae. 



3.2. Radial Velocities 



We computed the mean observed radial velocities of the ansae from 10 emission lines 
fitted with a Voigt profile to be -58.5 ±1.9 km.s"^ (E) and -47.9 ±1.6 km.s"^ (W). When 

using only the [Nil] and Ha lines (the HST images were taken in these lines) we obtain 
— 58.3±1.1 km.s^^ (E) and — 47.7±0.8 km.s^"*^ (W). The radial velocity differences between 
the E and W ansae are: Afja[Nii] = 10.6 ±1.4 km.s^-*^ and A^h = 10.6 ±2.5 km.s^^. We also 
separately computed the mean values of the velocity for the forbidden and for the permitted 
lines. They were also the same to within the errors. The observed radial velocities for 
each detected emission line are hsted in Table 3, and we also give the heliocentric and LSR 
velocities and their average values. 



3.3. Electron temperatures and densities 

The line intensities of [Nil] and [Sll] were measured by fitting Voigt profiles to derive 
the electron temperature and density in both ansae. The line intensities were corrected for 
reddening, assuming a recombination model B with T — 10^ K and rig = 10^ cm~^, and 
an extinction coefficient cp — 0.26 (Lame & Pogge 1996). The fine intensity ratios together 
with electron temperatures and densities derived from them are listed in Table 4. The latter 
results were obtained using the formulae from McCall (1984). The errors in the ratios were 
calculated from the S/N in each line. 
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Table 3: Radial velocities in km.s ^ for each emission line identified. 



eastern ansa 




western ansa 






jjine V ohs 


V hel 


^Isr 


^ obs 


^hel 




[iMij o/oo -oy.i/ 


-55.43 


-45.42 


-47.70 


-43.71 


-33.73 


Hel 5876 -56.42 


-52.68 


-42.67 


-48.51 


-44.56 




[01] 6364 -62.23 


-58.49 


-48.48 


-50.44 


-46.49 


-36.47 


[01] 6300 -61.19 


-57.45 


-47.44 


-49.05 


-45.10 


-35.08 


[Nil] 6548 -59.33 


-55.59 


-45.58 


-48.34 


-44.39 


-34.37 


Ha 6563 -56.69 


-52.95 


-42.94 


-48.46 


-44.51 


-34.49 


[Nil] 6583 -57.87 


-54.13 


-44.12 


-46.48 


-42.53 


-32.51 


[Sll] 6716 -57.68 


-53.94 


-43.93 


-45.34 


-41.39 


-31.37 


[SIX] 6731 -56.83 


-53.09 


-43.08 


-45.46 


-41.51 


Qi /in 
-oi.4y 


[AIII] 7134 -57.60 


-53.86 


-43.85 


-48.97 


-45.02 


-35.00 


Averages -58.5 ±1.9 


-54.76±1.9 -44.75 ±1.9 -47.9 ±1.6 -43.92 ±1.6 


-33.91 ±1.6 


Table 4: Derived line properties. 










Eastern ansa 












Ratio 


Observed 


Corrected 


Result 






[Nll](6548+6583)/5755 


111 ±15 


103 ± 15 


Te = 9000 ± 400 K 






[Sll] 6716/6730 


0.67 ± 0.04 


0.67 ± 0.04 


rie = 2600 ± 500 cm- 
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Western Ansa 












Ratio 


Observed 


Corrected 


Result 






[Nll](6548+6583)/5755 


120 ± 12 


110 ± 12 


Te = 8900 ± 400 K 






[SII] 6716/6730 


0.64 ±0.04 


0.64 ±0.04 


He = 1900 ± 300 cm- 
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4. 


Discussion 







Our radial velocities agree with those of Reay & Atherton (1985), who determined 
±6.2 km.s~^ using the [OI]6300 line. We get ±6.1 km.s~^ using the same line. 



-8- 



Liller (1965) computes IGmas-yr"^ for the angular expansion of the ansae; we obtain 
28 ± 8 mas.yr"^. This difference is probably due to the low resolution of the photographic 
plates, as Liller comments: "....the plate scale is such that photographic grain often competes 
with seeing as the limiting factor in image definition" . 

If the ansae were ejected at the same time from the CS and with equal & opposite 
velocities, we can determine the PM of the CS (relative to the point of ejection, assumed 
fixed in space) and the age of the ansae from our data. We find that the CS has moved 
0.96±0.04" since the ejection of the ansae, 920±260 yrs ago. This gives a proper motion 
for the CS of 1±0.5 mas.yr"^ or 5 mas between our two epochs, negligible cf. the motion of 
the ansae. 

In the gas of NCC7009, the sound speed is about lOkm.s"^, implying that the ansae 
move at supersonic speeds. 

The Te and Ug we obtain are similar to those found by others: Balick et al. (1994) quote 
a value of Tg = 8100K & Ug = 1000 cm~^; Bohigas, Lopez & Aguilar (1994) quote a value of 
Te = 9800K & Ue = 2300 cm'^, both for the W ansa. 

5. Conclusions 

We conclude that the ansae in NGC7009 -which is at a distance of 860± 340pc- move 
outward at the supersonic velocity of ± 112 km.s~^ near the plane of the sky (2.6 ± 1 degrees), 
have a dynamic age of 920±260 yrs, a proper motion of 28±8mas.yr~^, and that the CS 
likely moves at 1 ± 0.5 mas.yr"^. 
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